Pinning and Depinning of the Bragg Glass in a Point Disordered Model 
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The three-dimensional frustrated anisotropic XY model with point disorder is studied with both 
Monte Carlo simulations and resistively-shunted-j unction dynamics to model the dynamics of a 
type-II superconductor with quenched point pinning in a magnetic field and a weak applied current. 
Both the collective pinning and the depinning of the Bragg glass is examined. We find a critical 
current Ic that separates a creep region with unmeasurable low voltage from a region with a voltage 
V ~ I ~ Ic, and also identify the mechanism behind this behavior. It is further found to be possible 
to collapse the data obtained at a fixed disorder strength by plotting the voltage versus TI, where 
T is the temperature, though the reason for this behavior is unclear. 

PACS numbers: 74.25.Dw, 74.25.Qt, 74.25.Sv 



The behavior of elastic structures in the presence of 
point disorder is a profound question in condensed matter 
physics with relevance e.g. for charge density waves and 
vortex lattices. Due to the competition between the re- 
pulsive interactions, which favor a periodic structure, and 
both thermal fluctuations and quenched disorder, with 
the eS'cct to weaken this order, already the static prob- 
lem is a very difficult one. The common picture has for 
some time been that the quenched disorder turns the vor- 
tex lattice into a Bragg glass, which is a phase with alge- 
braically decaying correlationslli,i2i|, though some recent 
papers suggest a more complicated phase diagram[3iy|. 

The dynamical properties of a point disordered system 
in the presence of a strong driving force is an active field 
with several open questions and competing scenarios. In 
contrast, the behavior at weak fields, which is the subject 
of the present Letter, has been less debated, and the usual 
picture is that of a sharp depinning transition at zero 
temperature that is rounded through thermal activation 
at non-zero T . 

In this Letter we present results from dynamical simu- 
lations on a three-dimensional (3D) XY model that give 
evidence for a different and more interesting behavior. 
We find a critical current Ic that separates a creep re- 
gion with unmeasurably small voltage from a region with 
a rectilinear behavior. This finding is shown to be re- 
lated to a memory effect in the system and is thus man- 
ifestly an effect of the system being out of equilibrium. 
The obtained voltage may furthermore be collapsed in an 
unexpected way. Our simulations are in many respects 
similar to Ref . |5|, but our longer simulation times make it 
possible to probe the behavior with considerably higher 
precision. 

The Hamiltonian of our 3D XY model is 
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links between nearest neighbors. The vector potential 
Ai^ is chosen such that V x Ai^ = 27r/z. / = 1/45 
gives one vortex per 45 plaquettes, which is less dense 
than in earlier simulations and is chosen with the aim to 
reduce the effects of the discretization by the numerical 
grid. The system size is 45 x 45 x 32 which gives 45 vor- 
tex lines. The anisotropy is chosen as Jiz/J = 1/40 and 
the disorder is introduced as point defects as in Ref. [tI 
through a low density {— 1/180) of plaquettes in the x-y 
plane with four weak links. The weak links have cou- 
pling (1 — p)J instead of J. In this study of the vortex 
solid we use disorder strenghts up to p = 0.34. By mea- 
suring the energy of the four links around a vortex both 
for defect positions and ordinary plaquettes we find that 
the pinning energy — the energy needed to displace a vor- 
tex from a defect position to an ordinary plaquette — is 

To check for collective pinning we used standard 
Metropolis Monte Carlo simulations with fluctuating 
twists 6] ; the driven vortex solid was studied with RS J 
dynamics with fluctuating twist boundary conditions _8| . 
The equations of motion are as in |9j and to integrate 
them we used a second order Runge-Kutta method with 
a dimensionless time step of At — 0.1 and typically (1- 
4)xl0^ units of time per run. The twist variable in the 
z direction is kept fixed, A^ = 0. 

Since the motion of the vortices is linked to the change 
of the twist variables, our main quantities are extracted 
from the behavior of the twist variables A^^ and Ay. In 
the absence of an applied current A^ will perform a ran- 
dom walk and the linear resistance is directly propor- 
tional to the diffusion constant, 



Ri^, ^ ^-^{[A,{t) - A,{0)f 
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Here 0i are the phase angles, the A^ are twist variables 
in the three cartesian directions, and the sum is over all 



A vanishingly small R means that the vortex lattice es- 
sentially stays flxed. An applied current gives a force on 
the vortices, and their motion means a change in the twist 
variables, which is seen as a voltage across the system. 



The voltage per link is given by 



Lp dt 



To monitor the order we measure the structure factor 
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where V = L^LyLz, and define Si as the algebraic aver- 
age of S'(k) for the six smallest reciprocal lattice vectors, 
k 7^ 0. We also trace out the vortex lines; the notation 
Yi{z) is used for the position in the x-y plane of vortex line 
i at plane z. The mean-squared fluctuation of the vortex 
lines (w^), is also determined. The presented data are 
for a single disorder realization, but similar results have 
been obtained with other disorder realizations. 

We now turn to the equilibrium dynamics of the point 
disordered system and show the temperature-dependence 
of the linear resistance Ri together with the structure 
factor Si in Fig. ^ The onset of Ri coincides here with 
the vanishing of the structure factor and the direct in- 
terpretation is that the vortex solid is collectively pinned 
by the disorder whereas the vortex line liquid is mobile. 
The data in Fig. Q] is for rather strong disorder, p > 0.2. 
At weaker disorder, p = 0.1 (not shown), the VL for the 
size 45 X 45 X 32 is floating whereas collective pinning is 
restored for a bigger size, 90 x 90 x 32. Based on the 
data in Fig. ^ together with some additional simulations 
at larger p, we obtain the phase diagram in Fig. El The 
solid symbols are the parameter values where the dy- 
namic simulations described below have been performed. 
The dashed region is where we expect effects of the nu- 
merical grid to be important, as discussed further below. 

After having established the existence of collective pin- 
ning we turn to our main subject, which is the dynam- 
ics of a weakly driven system. At zero temperature the 
vortex lattice depins from the numerical grid at the crit- 
ical current density Jq « 0.13, but our focus is on de- 
pinning from the disorder potential which takes place 
at current densities more than two orders of magnitude 
smaller than Iq ■ In the simulations the current is applied 
in the diagonal direction with Ix = ly = I/V^ since 
this gives the force // along — v = (x — y)/-\/2 which is 
along one of the symmetry directions of the vortex lat- 
tice, see inset of Fig.QIa). Likewise the voltage per link 

is y = (14 + Vy)/V2. 

Figure 13 shows the I-V characteristics at T = 0.24 
with four different pinning strengthsjlOJ. In the pure 
system, p = 0, the voltage is proportional to the current, 
V = RqI. The main effect of non-zero pinning strength 
is to shift the voltage down towards zero. As this shift 
downwards only depends weakly on /, the slope is almost 
independent of p. In this way we get a critical current 
density Ic that separates the creep region with an unmea- 
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FIG. 1: Collective pinning in the vortex solid. Panel (a) 
shows the linear resistance Ri = |(-R;,a; -I- Ri,y) determined 
from Eq. Q with t = 16384 from data obtained with MC 
simulations. (The non-zero Ri at low T is an artifact of a 
rather small t together with the fluctuations of A^ around a 
fixed values.) The inset is a top view of a configuration of vor- 
tices; the darker points indicates a higher density. Panel (b) 
shows the structure factor Si ; the inset is an intensity plot of 
ln5'(k). Both insets are for p = 0.30 and T = 0.30. Note that 
the resistance, and thereby the mobility of the vortex solid, 
decreases rapidly as the temperature drops below melting. 



surably low voltage at /< Jc from the moving region at 
I > Ic with V r^ I - Ic- 

This behavior suggests the existence of a force which, 
much like a friction force, has to be overcome to set the 
vortices in motion and also reduces the velocity of the 
moving system. For a system in equilibrium one would 
expect the random pinning to contribute with a number 
of random forces which tend to cancel one another out, 
but the slow drift of the vortices opens up for other possi- 
bilities as e.g. memory effects. The idea behind our sug- 
gested mechanism is that some of the defect vortices ini- 
tially stay fixed when the unpinned vortices slowly move 
in the direction of the force. This leads to distortions 
of the elastic vortex lines (see inset of Fig. OJa)) and a 
gradual build-up of a force /d on the defect vortices. This 
also gives a reaction force fd from the defect vortices on 
the vortex lattice; the motion of the vortex lattice should 
then be proportional to fi — fd- 

A basic ingredient in this scenario is that the vortices 




FIG. 2: Phase diagram for the simulated model. The open cir- 
cles connected by a solid line is from determinations of melting 
temperatures from the vanishing of the structure factor; the 
continued dashed line is a sketch of a possible behavior at low 
temperatures. The "Floating solid" region at low p is where 
the disorder strength is too weak to pin the vortex solid. This 
is a finite size effect; pinning is restored at a larger size. The 
filled circles show the location in the phase diagram of the 
runs in Fig. Olat T = 0.24 and the runs in Fig. glat p = 0.32. 
The dashed region at low temperatures is where the numerical 
grid is expected to affect the results. 



at the defects on the average lag after the unpinned vor- 
tices. To test this idea we introduce a way to measure 
the average displacement of the defect vortices. To dis- 
tinguish between ordinary vortices and defect vortices (or 
rather defect positions) we introduce the defect operator 
D(r, z) which is unity at a defect and zero otherwise. 
The average position of all non-defect vortices belonging 
to the same vortex line may then be written, 



r, -[1-D(r,(z),z)]r,(z), 



and as a measure of the lagging behind of the defect 
vortices, we introduce the displacement as the sum of 
their deviations from the average vortex line positions, 






(r,(z) -r,) D{y,{z),z) 



(3) 



Note that this is only the displacements relative to the 
vortex lines. An estimate of a total displacement should 
also include effects due to a vortex line or a set of vor- 
tex lines lagging behind the rest of the vortex lattice. 
This is however considerably more difficult to estimate, 
and as we will see shortly it seems that the expression 
above contains the dominant contribution. Nevertheless, 
d is only an approximate estimate — most likely a lower 
bound — of the total displacement. 

We now expect the displacement vector to be non-zero 
and point in the direction opposite to the force. Fig- 
ure El^b) which shows how d • v depends on the applied 
current gives strong evidence for the suggested behavior; 
the defect vortices on the average lag after the ordinary 
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FIG. 3: Results at T = 0.24 and several different p. Panel (a) 
shows that the effect of a finite p on the I-V characteristics 
is to shift the data down towards zero. This gives a creep 
region with very low voltage at low currents, I < Ic, and a 
region with a rectilinear behavior at larger /. The inset which 
is a vortex line extending along z and moving in the direction 
of the applied force, illustrates the suggested mechanism: the 
vortices at defect positions on the average lag after the moving 
vortex lattice. Panel (b) gives the displacement, c.f. Eq. Q in 
the direction opposite to //. Panel (c) which is a calculation 
of the voltage from Ro of the pure system and Id from the 
measured d and Eq. , shows that the suggested mechanism 
captures the essential elements of the dynamics. 



vortices. For a qualitative comparison one would like to 
compare the magnitude of this defect force fd with // 
from the applied current. We do a similar but more di- 
rect comparison by instead estimating Id, the current in 
the system due to the displaced vortices. Recall that the 
creation of a vortex pair in a 2D model with a separation 
of d lattice constants in the y direction creates d rows 
in the x direction where the phase angle rotates by 27r. 
This gives a total current /tot ~ dsin(27r/La;) « 2nd/ L^- 
Taking this over to our 3D system we conclude that a 



displacement d corresponds to a current density 
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At I < Ic we expect Id = I and it is then possible to ob- 
tain the displacement d from Eq. 10}. This is the dashed 
line in Fig. I^Jb). The data for I < Ic agrees well with 
this line. For the behavior at higher currents we plot in 
Fig.EIc) Vd = Ro{I—Id), which is an estimate of the volt- 
age based on the measured displacement d, together with 
Rq from the pure system. The agreement with Fig. ISfa) 
is good and shows convincingly that our simple descrip- 
tion in terms of the displacement of the defect vortices 
captures the essential elements of the dynamics. 

We now fix the pinning strength, p = 0.32, and exam- 
ine the I-V characteristics at several different tempera- 
tures from T = 0.08 to T = 0.30. The results are shown 
in Fig. 0Ja) and give evidence for a critical current that 
increases with decreasing temperature. In a further anal- 
ysis we have found that it is possible to collapse the data 
by plotting the voltage versus the combination T/, see 
Fig. Hb). This gives Ic oc 1/T and ViTJ) = xil/Ic) 
where x is a scaling function. To check if this behavior 
is general or depends crucially on our way of introduc- 
ing the disorder as a low density of very strong defects, 
we have done additional simulations with random cou- 
plings as in Ref. [llj with p = 0.08. These simulations 
give support for the same kind of scaling collapse. 

At low temperatures one expects an effect of the dis- 
crete lattice would be to slow down the dynamics. This 
is precisely what is seen at T < 0.12 (the data fall below 
the common curve) and this data is therefore discarded 
to produce the collapse in Fig. ^b) . To study the effect 
of the discretization we return to the behavior of the lin- 
ear resistance in a clean system. At high temperatures 
the linear resistance decreases slowly with temperature 
down to T « 0.22 where it starts falling rather rapidly 
down to i?; = at T = 0.14. This shows that the effects 
of the discrete grid are significant up to T « 0.22. Since 
the discretization is expected to be unimportant if the 
typical wandering of the vortex line is much larger than 
the lattice spacing, we examine the mean-squared fluctu- 
ation which is (u^) = 1.03 at T = 0.22. It also turns out 
that the data at low T that were discarded in Fig. 2{b) 
all have (^m^) < 1 whereas (w^) > 1 at the temperatures 
included in the figure. This suggests (w^) > 1 as a crite- 
rion for significant effects due to the discretization (the 
dashed region in Fig. ^ , and gives support to our claim 
that the failure to collapse the data at low T actually is 
an effect of the discretization of the underlying grid. 

The origin of the collapse in Fig. Efb) is surprising as 
several quantities that could be expected to affect the 
mobility change considerably across this temperature in- 
terval. Examples are the fraction of vortices on defect 
positions that decreases from 0.078 at T = 0.14 to 0.054 
at T = 0.30, the structure factor that decreases from 0.5 
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FIG. 4; The voltage at p = 0.32 and several different tem- 
peratures. Panel (a) gives evidence for a critical current that 
increasese with decreasing temperature. Panel (b) shows that 
the voltage may be collapsed to an impressing accuracy by 
plotting the data versus TI. 



to 0.3, and the resistance in the absence of disorder which 
behaves as Ro{T) ~ T^/^. 

To conclude, we have obtained the I-V characteristics 
of a point disordered frustrated 3D XY model and found 
that the effect of disorder is to give a critical current that 
separates the creep region from a region with V ^ I — Ic- 
The reason for this behavior is a pinning force which 
appears to be due to a memory effect; a fraction of the 
defect vortices lag after the rest of the vortex lattice. 
Finally, it is shown that the voltage may be collapsed by 
plotting the voltage versus TI, though the origin of this 
behavior as yet not is fully understood. 
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